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ABSTRACT 

The universe may have undergone a superfluid-like phase during its evolu- 
tion, resulting from the injection of nontopological charge into the spontaneously 
broken vacuum. In the presence of vortices this charge is identified with angu- 
lar momentum. This leads to turbulent domains on the scale of the correlation 
length. By restoring the symmetry at low temperatures the vortices dissociate 
and push the charges to the boundaries of these domains. If we-(phenomenolog- 

T t* -ft* * ' 4 ~~ ' ** ^ ' 

ically) scale our model to very low energies, wecan incorporate it in a late-time 
phase transition and form large scale structure in the boundary layers of the 
correlation volumes. The novel feature of our model lies in the fact that the 
dark matter is endowed with coherent motion. We^laborate on the possibilities 
of identifying this flow around superfluid vortices with the observed large-scale 
bulk motion. If this identification is possible, then we- can make the definite 
prediction that a more extended map of peculiar velocities would have to reveal 
large-scale circulations in the flow pattern. 
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1. Introduction 

In recent years it has become increasingly evident that the universe may 
be structured on very large scales, extending to distances of 100 h~ l Mpc and 
possibly beyond. A variety of different theoretical models have been proposed 
to explain these observations. They include the ‘standard’ scenario [1] (hot or 
cold dark matter with a fluctuation spectrum from inflation), the ‘seed induced’ 
models (cosmic strings [2], textures [3], etc.), and late-time (lifter decoupling) 
phase-transitions in the cosmic evolution [4-7]. Although these models differ 
from each other in many respects, and have their own difficulties, most of them 
share a common point: they try to explain the observed large-scale peculiar 
velocities by gravitational infall [8]. This assumption is usually thought to be 
reasonable, as every peculiar velocity, not supported by a gravitational perturba- 
tion, is redshifted during expansion [9]. In fact, this consideration is the backbone 
of one of the most powerful tools in tracing the dark matter in the universe [10]. 
It directly links the observed peculiar velocities with the underlying large-scale 
density perturbations. Based on this very connection, there have been indications 
for some time, that ft (the ratio between the total density to the critical density) 
may be ^ 0.5 [10,11]. A significant improvement has recently been made in the 
statistical analysis of the observational data and it is now claimed that ft <: 0.5 to 
a 95 % confidence level [12]. As primordial nucleosynthesis constrains the maxi- 
mum amount of baryonic matter in the universe to ftbaryons — (0.011 — 0.048)h" 2 
[13], this result leads to the first strong observational evidence for the existence of 
non-baryonic dark matter. There clearly still exist observational and theoretical 
uncertainties [14] in the above cited results, and it may well be that the constraint 
^ eventually ‘softens’ somewhat. In the meantime we are, however, called upon to 
look for alternatives, which could circumvent such a conclusion. 

A model for non-gravitationally induced velocities was already proposed a 
few decades ago [15] in an attempt to explain the angular momentum of galaxies 
and the origin of magnetic fields. It was realized that the cosmic plasma may 
be in a state of turbulence during the radiation dominated era and thus lead to 
vorticity in the flow. This scenario is, however, ruled out as it leads to structure 
formation too early on [16]. 

We can, however, evade this problem by letting the universe undergo a 
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superfluid-like phase in its late evolution (at a redshift of z = 10 — 10 3 ). Po- 
tential flow around superfluid vortices could then directly evolve into coherent 
bulk motion, at the same time as the formation of structure occurs. All the neces- 
sary ingredients for such a scenario may be found in a model, recently proposed 
by Dodelson and Widrow [17] to account for the baryon number in a baryon- 
symmetric universe. In this respect, our scenario may be viewed of originating 
to some extent in a high energy theory. However, in order to make contact with 
the actual cosmological data, we have to rescale the model to very low energies 
and to require tiny coupling constants and mass parameters and therefore end 
up with very light particles. This is a well known drawback for most late-time 
phase transitions, though it has been solved in some cases, as for example in the 
schizon model [18]. Consequently, with respect to fixing parameters, our theory 
is clearly phenomenological [19]. 

This poses, however, a basic problem. Our goal is to find a way to generate 
large-scale velocities and, if possible, to challenge the conclusion about the exis- 
tence of non-baryonic dark matter, i.e., is non-baryonic dark matter necessary? 
However, in our ‘low-scale’ model we eventually end up with very light particles, 
which could not possibly be identified with baryons. We therefore again introduce 
non-baryonic matter. Our only achievement would then lie in our mechanism of 
producing coherent motion in the dark matter. This then still leaves us with 
the difficulty of transferring the vortical motion efficiently into the baryonic com- 
ponent. Alternatively, however, we may view our model as a first attempt of 
introducing a late-time baryonic superfluid phase into cosmology. This paper 
then presents an outline of the maun ingredients and ideas necessary for such a 
scenario. It is then hoped that a realistic model can eventually be found in a 
more comprehensive context. 

Assuming our scenario to be correct, we make a definite prediction for future 
observations: a more extended map of peculiar velocities would have to reveal 
large angular motions (circulation), with the currently observed radial velocities 
being only a part of the more general picture. 

The paper is organized as follows. In the next section we briefly review the 
idea of injecting nontopological charge into a coherent background. In sect. 3 
we introduce ‘vortices into the system and show how the scenario changes signifi- 
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cantly. We then try to make contact with observations. Sect. 4 contains a short 
summary. 

In the final stage of this work we learned about a preprint by Davis [20], 
in which he also addresses the idea of generating large-scale structure out of a 
coherent background and a network of defects. His scenario and results differ, 
however, from ours. Nevertheless, in order to accomplish enough charge density, 
the model of ref. [20] suffers from essentially the same ‘small-parameter’ problems. 

2. Coherent State 

Our effective lagrangian incorporates a complex scalar field (f> with a U(l) 
symmetry. (For the moment we will closely follow the presentation by Benson 
and Widrow [21].) The temperature dependent potential of the theory reads 

V(*,T) = (!) 

where m,M are mass parameters and A is a dimensionless coupling constant. 
This potential leads to spontaneous symmetry breaking at high temperatures and 
to a restored symmetric phase at low temperatures. It is assumed that in the 
broken phase an external field (‘external’ only in the framework of the effective 
theory) decays into U(l) charged particles. As the theory is U(l) symmetric, the 
produced charge must be compensated by the same amount of ‘anti-charge’. The 
charge (or anti-charge) may however be hidden in the vacuum by way of putting 
the background into a coherent state. This is similar to the interior of a Q-ball 
[22]. (In the original model of ref. [17] baryon number was identified with the 
U(l) charge. The external field then decayed into free baryons and a coherent 
vacuum, which carried the anti-baryonic charge.) 

Let us first assume a spatially homogeneous field, 

m = e ‘ e{,) ■ < 2 ) 

The conserved Noether current, j M leads, for the coherent 
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state (2), to the charge density 


n 4> = j° = (3) 

with (3 = 0. 

The equations of motion in a Friedmann-Robertson- Walker universe with 
scale factor a read 


a 8V • , 

v + 3 — v + — - v 0 2 = 0 , (4) 

a ov ' 

® + ( 3 f + 2 z) ® - 0 • < 5 > 

Integrating eq. (5) gives 

n <t> = v 2 0 = VT* , ( 6 ) 

where we have assumed an adiabatic expansion, a <x T -1 ; here q is a dimension- 
less constant. The total charge inside a comoving volume, Q = u 2 (3V, is therefore 
conserved [23]. 

Let us first review the main points of the original scenario [17,21]: At high 
temperatures the vacuum is in its broken state, u ^ 0. As the external field decays 
and produces charged particles, the charge of the anti-particles (henceforth called 
^-particles or just particles) is absorbed by the vacuum in form of the coherent 
state (2) with (3 ^ 0; the Goldstone boson starts ‘rotating’ around its vacuum 
manifold. As the universe expands and cools, the background undergoes a first- 
order phase transition to its symmetric state at T c cs ^/3/A m. However, a 
symmetric vacuum with u = 0 cannot support any charge and the coherent 
background field must either decay into free particles or, alternatively, stay in 
a broken phase, an effect well-known from ‘charged-induced’ phase transitions 
[24]. It is expected that in the latter case only localized regions will be trapped 
in the broken state and carry all the charge and therefore show up as isolated 
nontopological solitons [21]. 
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3. Superfluid Vortices 

In the scenario outlined above we have neglected the possibility that the sys- 
tem, at high temperatures, could include topological defects. Based on causality 
arguments there is, in fact, no reason (at least in the absence of inflation) to rule 
out global strings. The presence of strings, however, changes the above scenario 
significantly. Around a globed string we can no more describe the ‘charged’ vac- 
uum by the spatially homogeneous field (2). Instead, in the presence of a straight 
static string along the z-axis, <f> must be modified to 


ff, = e »("*+©(0) ; 

\/2 


( 7 ) 


r and 6 are polar coordinates in the x-y plane and n is the winding number of 
the string, i.e. the topological charge. If we neglect the radial dependence of 
v (an approximation which is justified outside the string’s core), then the total 
nontopological charge (for a given /?) in a comoving volume around the string is 
unchanged. However, the energy-momentum tensor, diagonal for the string-free 
background (2), develops now the non-diagonal term Tg and therefore angular 
momentum. The axisymmetric system admits the Killing vector (Jg)*' and hence 
the conserved current 



The magnitude of the angular momentum (in the z-direction), for our coherent 
state (7), becomes 



( 9 ) 


J is conserved in a finite volume as long as no angular momentum flows through 
the boundaries. According to (9), angular momentum is directly related to the 
nontopological charge in the system. (In the following we set n = 1.) The 
angular momentum is correlated over a volume = 7 r£ 3 , where £ is of order 
the interstring distance. (We picture V( as a cylinder of radius £ around a string 
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segment of length £.) Injecting charge into the vacuum in the presence of strings is 
therefore equivalent of endowing the correlation volume with angular momentum. 
If £ <C H~ l we end up with ‘turbulent patches’ on scales £. 

During the decay of the external field angular momentum is of course con- 
served, thus J of the coherent field is exactly canceled by the angular momentum 
of the free ‘anti-^-particles’. However, due to the lack of any coherent motion 
of these particles, they freely stream into neighbouring correlation volumes and 
thereby smear out any initial angular momentum in their component. 

The actual size of the correlation volume is fixed by the dynamics of the 
strings. Before injecting the charge into the vacuum at a temperature T in, the 
global strings move relativistically and, once inside the horizon, decay due to 
Goldstone boson radiation in a few oscillation times [25]. The correlation length 
of the string network will therefore scale with the horizon, £ ~ t. At T < T{ n 
the strings are imbedded in a background. Depending on the actual charge 
density, the string motion may become highly damped - the global strings es- 
sentially behaving like superfiuid vortices [26] (or spinning strings [27]). Let us 
be more specific: the time-dependence of the Goldstone mode in eq. (7) leads 
to a Lorentz-like force (the Magnus force in the terminology of fluid dynamics 
[28]), counterbalancing the string’s tension. The tension of a global string is 
given by /tension — {™ 2 /O ln (Z/ 6 )i where 8 is the core radius. The Magnus force 
reads [26,27] /Magnus = 7 ti/ 2 / 37 v; t; is the string velocity and 7 = (1 — -u 2 ) -1 / 2 . 
Neglecting the radiation field and the ‘induced Magnus force’ [26], we can set 
/tension = /Magnus and reduce for the string’s velocity 

1 

( 10 ) 


V — 


1 + 


($)■ - 


MS) 


with Tf = (3 1 . Clearly, if /3 — ► 0 (no charge) then the global strings move highly 
relativistically. However, if 

- » HW . (u) 

T f 

then v ~ By taking the limit Tf —* 6 we recover the nonrelativistic 

behaviour of vortices in superfiuid 4 He [28]. 


7 


Whether (11) is realized depends on rj and 6 and hence on the explicit 
solution to the equations of motion. For the moment let us assume that eq. (11) 
applies for temperatures below T sv , T sv & Tin- For T < T 3V the globed strings 
mimic superfluid vortices and are essentially frozen into the background. The 
string network is then only conformally stretched, £(<) = a(t)£ sv (£ S v = 

During the subsequent expansion of the universe, string length accumulates inside 
the horizon and the Hubble volume becomes filled up with an increasing number 
of correlation volumes [29]. 

In order to check whether (11) can indeed be satisfied we have to compute 
Tf and 6. Substituting (6) into eq. (4) allows us to rewrite the v equation in the 
form 


+ 3 — if + 
a 


dV e 


e f f 


dv 


( 12 ) 


with 

V„,{v,T) = V>,T) + ^ • (13) 

Instead of solving eq. (12) explicitly, we follow Benson and Widrow [21] and 
assume that u traces the minimum of V e ff- We therefore look for a solution to 
dV e ff/d v = 0, i.e. 

y 3 + 3p y + 2q = 0 , (14) 


where y = (£) 2 , p = ^ and q = §p - {jfcsV 2 + =£)/2 tj 2 . By defining u± = 
[— q ± (g 2 +P 3 ) 1 / 2 ] 1 / 3 the solution to eq. (14) is given by y = u+ + «-, which can 
then be inverted to yield u = i/(T). The effective potential of eq. (1) has been 
derived using a perturbation expansion [21] and requiring (^r) (773-) = A 2 . 

We may therefore neglect p and find 

/ T \ 3*/4 

u(T) ~ 2( 2 ~ fc >/ 4 u ( , (15) 

where v = (m 2 il/ 2 /4A) 3 / 4 , T = (m.*M 2 / 4A77) 1 / 6 and k = 1 or 2 for T T and 
T T, respectively. Inserting v from eq. (15) into eq. (6) yields 


T f = / r 1 ~ ( 2 ( 2 -*)/ 2 V 2 h T ^ 2 ) T W 2)-3 . ( 16 ) 


Note that in a radiation dominated universe, rf grows more slowly than the 
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horizon, whereas in a matter dominated universe (the case we are interested in), 
rf ~ H~ l in the high temperature phase. But even during matter domination, 
eventuaUy £/ rf 1. (The fastest scaling of £ is ~ t, the slowest ~ a.) This 
means that there exists a temperature T sv , below which eq. (11) is satisfied and 
the strings become ‘slow-movers’. (For completeness we could easily estimate 
the core radius by £ ~ m^ 1 , where is the mass of the excitations around 
the spontaneously broken u.) Our only worry is that T sv > T c . Whether this is 
possible clearly depends on the parameters of the theory. 

The process of string accumulation will continue until the background un- 
dergoes the first-order restoring phase transition at a temperature T c and the 
vacuum is driven to its symmetric state. Now we have to recall that the cen- 
tred core axis of the vortices is already in the true vacuum. During symmetry 
restoration the vortices are acting as nucleation sites [30], Their cores expand and 
thereby convert false to true vacuum. However, as the symmetric background 
cannot support any nontopological charge, the charge is pushed out to the bound- 
aries of V^. This will continue until an equilibrium configuration is reached, in 
which the pressure from the particles (charge) balances the pressure from the 
potential difference. Let us assume that this process is short compared to H~ l 
and therefore neglect the expansion. Consider a cylindrical correlation volume 
of height £ and base n£ 2 around the straight and static vortex. After symmetry 
restoration the vortex core expands to a radius r c — £ (0 < a < 1) and 

pushes out all the charge to the cylindrical layer between r c and £. Neglecting 
the gradient energy of i/.at r e (where v changes from 0 to v ^ 0) the total energy 
inside is 

E = 2*( jrdr ( \ fp + V(„) + \ £ ) . (17) 

Tc 


The first term results from the kinetic term of the Goldstone boson, V{y) is the 
potential energy in the outer layer (inside the expanded core V(0) = 0) and the 
last term yields the usual logarithmically divergent energy per unit string length 
of a global string. Since T < T sy and hence (3 £ -1 , we can neglect the last 
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term in (17) to obtain 


E = ( \ v 2 P 2 + V(v) ) S( , (18) 

with 5 = Q7r^ 2 . In order to evaluate r c we have to minimize the energy in a 
specific Q-sector, i.e. for given nontopological charge. Substituting /? = Q/t/ 2 S£ 
into eq. (18) and minimizing with respect to S yields 

a ~ ( 2 v 2 V(u) ) (19) 

and the energy 



Here is the charge density in the absence of vortex dissociation. Eq. (20) is 
reminiscent of a nontopological soliton [22], In fact, it is precisely what we are 
looking for, namely a soliton-like configuration in the outer layer of the cylindrical 
correlation volume. The energy in eq. (20) must still be minimized with respect 
to u in order to find the field value in the soliton layer. The result is z / m j n = 
( y/Z/2)M . Inserting t /min into eqs. (19) and (20) then yields r c and the energy of 
the soliton layer in terms of the fundamental parameters and the temperature. 
We would like to have a — ♦ 0 and thus r c — ► f in order to push the charge as far 
out as possible. 

Though we have been looking for an equilibrium soliton-like configuration in 
the outer layer of V£, there is in fact no reason to believe that such a state will 
ever exist. During the process of vortex dissociation the coherent field becomes 
excited and may easily decay into free particles, fragment into spherical Q-balls, 
or a mixture of both. This does not contradict with topological arguments, as v 
can vanish along a whole radial direction and thus allow the coherent state with 
nonzero winding number to ‘unwind’. Our only assumption is that the charge is 
pushed out significantly before the fragmentation or decay process begins. If r e 
is of the order of £ then we end up with high charge concentrations in the layers 
of radius ~ £ and charge-free space in the interior of the expanded vortex cores. 
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Identifying the charge (the (f) particles) with the daxk matter in the universe, 
this scenario leads naturally to structure formation in the high density regions. 
The (curved) sheet-like structures of correlation length £ are then intercepted by 
voids of volume £ 3 . The novel feature of this scenario resides in the fact that no 
matter what form the charge ends up after the decay of the coherent cylindrical 
configuration, angular momentum must be conserved. This naturally provides 
the dark matter with coherent angular motion. 

So far our scenario fits nicely into the model of refs. [17,21]. We now make 
contact with observations and try to identify the large-scale streaming velocities 
of v ~ 600 km/sec on a scale of £o ~ 60 Mpc with the coherent angular motion 
around the dissolved vortex lines. If most matter exists today in form of free 
<t> particles, then their number density is given by = p/m. By equating the 
charge, Q = with the magnitude of the angular momentum, J ~ M. £o v ~ 

P v, we obtain 

m ~ ~ 5- 10 -29 eV . (21) 

So v 

We need a very light and extremely abundant (n^ ~ 10 32 cm -3 ) particle in or- 
der to account for the observed bulk motion. Note that the ‘specific angular 
momentum’, £v, is constant during the expansion. 

For T < Tg V the correlation length is conformally stretched. Requiring £o ~ 
60 - 90 Mpc fixes the onset of the superfluid epoch to a redshift of z sv ~ 500 
- 1000. This, together with T e < T sv and eq. (21), fixes A and M. As 

expected, the parameters tire small. However, as we have already emphasized in 
the introduction, we view our model as a phenomenological one, having the same 
difficulties as most late- time phase transition models [31]. 

How large are the induced inhomogeneities in the microwave background 
radiation? The largest $£ is probably generated by the high density layers just 
after their formation, thus at a redshift ~ of the dissociation process. A 
rough upper limit can be borrowed from the soft domain wall scenario [6]: 

( 22 ) 

where l is the layer thickness and p the matter density inside the layer. Tak- 


„ ,2 
~ G P 1 , 
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ing p = p{i/l) 2 ( p denotes the appropriate density without vortex dissociation, 
homogenous throughout space) we estimate 

ST 3 

"jT — ( Bob ) ( 1 + ^dis ) — 5 * 10 5 ^layex ( 1 + z dis ) > (23) 

where fli ayer = p/ pc and p c is the closure density. The Hubble constant is taken 
to be Ho = 100 km/sec/Mpc. As the density perturbations are already nonlinear 
after dissociation, we can easily take Zdis to be as small as 5 - 10. Furthermore, 
there is no reason why a critical density should be ‘squeezed’ into the soliton-like 
layer. On the contrary: the total mass density can be below the critical density, 
say fitotal ~ 0.2, and fli ayer only a small fraction of this total density. This would 
lead, according to (23), to distortions in the CMBR below observational limits. 

4. Conclusions 

We have presented a scenario in which the universe experiences a superfluid- 
like phase during its late evolution and ends up with large-scale streaming ve- 
locities in the dark matter component. This coherent motion is not supported 
by any gravitational perturbation, but is a direct consequence of the potential 
flow around superfluid vortices. This is especially appealing in light of the recent 
analyses of gravitationally induced peculiar velocities, which claim evidence for 
non-baryonic dark matter. If we could find a more realistic theory for baryonic 
superfluidity, then we may be restrained of making such a firm prediction about 
the existence of non-baryonic dark matter. Furthermore, such a superfluid epoch 
would have a strong imprint on the velocity map and should reveal angular mo- 
tions (circulations) on very luge scales, beyond the currently probed distances. 
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